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Abstract 



Intermediate mass fragment (IMF) formation in the 12 GeV proton induced reaction 
on Au target is analyzed by using a combined framework of a transport model 
(JAM/MF) and a newly developed non-equilibrium percolation (NEP) model. In 
this model, we can well reproduce the mass distribution of fragments. In addition, 
the sideward peaked angular distribution would emerge under the condition that 
the fragment formation time is very short, around 20 fm/c. Within this short time 
period, the un-heated part of the residual nucleus is kept to have doughnut shape, 
then the Coulomb repulsion from this shape strengthens the sideward peak of IMF. 
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1 Introduction 



Multifragment formation from excited nuclei has attracted attention in this 
decade. It is closely related to the properties of nuclear matter, especially to its 
liquid-gas phase transition. For example, it has been proposed that multifrag- 
mcntation would be caused by the volume instability of nuclear matter in the 
spinodal region [1,2], in which nuclear liquid and gas phase rapidly co-exist. 
The above idea on the fragment formation has received renewal interest, since 
recent high-quality heavy-ion experiments provide us of strong evidences of 
the first-order nature of the nuclear liquid-gas phase transition [3-5]. It is a 
big challenge for theorists to elucidate the relation between the fragment for- 
mation and the equilibrium properties of nuclear matter, including its phase 
transitions. In heavy-ion collisions, however, non-trivial roles of strong col- 
lective flow make it difficult to relate fragment formations and the nuclear 
matter properties directly. Thus it is also necessary to investigate fragment 
formation phenomena in light-ion induced reactions, in which the flow effects 
on the fragmentation are expected to be smaller. 

It has been known that the fragment mass distribution in high-energy proton 
induced reactions follows the power law, which signals the statistical formation 
of fragments in a critical but equilibrium condition. Energy spectra of IMFs 
also support the statistical fragment formation picture. In a recent KEK ex- 
periment of 12 GeV proton induced reaction on Au target, the energy spectra 
of IMFs such as O, F, Ne, and Na are measured, and it is found to be well 
fitted with a thermal model with Coulomb barrier. 

Contrary to the mass and energy distribution of fragments, the angular distri- 
bution of IMFs contradicts to a picture of the statistical fragmentation from 
an equilibrated system. It is also long known that the angular distribution 
of Intermediate Mass Fragments (IMFs) becomes sideward peaked at proton 
incident energies above Ep = 10 GcV [6]. In the above KEK experiments, the 
angular distribution of IMFs such as O, F, Nc, and Na are shown to have a 
sideward peak around 70 degrees in the laboratory frame. The same exper- 
imental conclusion is also obtained in a BNL-AGS experiment. In the AGS 
p-|-Au experiment [7], the laboratory angular distribution of C isotopes be- 
comes to have sideward peak as the incident energy of proton go beyond 10 
GeV. This sideward peaked feature has been considered to be a mystery; it 
contradicts to a naive picture of statistical fragment formation from a spherical 
equilibrated residual nuclei. 

In a statistical picture of fragment formation, the angular distribution of frag- 
ments is generally forward peaked. At low energies, the incident proton is 
absorbed into the target, and forms a thermalized compound nucleus having 
the incident proton momentum. At high energies, the incident proton collides 
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with nucleons inside target nucleus (cascade process) in the first stage of reac- 
tions. These successive two-body scatterings give momentum transfers to the 
nucleons in target, and some of them are emitted from the residual nucleus. 
The momentum transfer to the residual nucleus is in the forward direction in 
average. Thus in both cases, the residual excited nucleus should have a forward 
momentum in average. If this residual compound nucleus decays statistically, 
produced fragments should have forward peaked angular distribution. 

It is expected that these fragment angular distribution anomaly would be due 
to a new fragment formation mechanism, which becomes evident at higher 
incident energies. This new mechanism would be related to a new property 
of excited nuclei produced in GeV proton induced reactions. For the under- 
standing of this sideward emission of IMFs, several characteristic mechanisms 
have been proposed so far; the shock wave propagation [8] , the strange shaped 
nucleus formation [9,10] which is also predicted to be formed in heavy ion 
collisions [11,12], and recoil effects of particle scatterings [7]. However, there 
is neither decisive conclusion nor satisfactory explanation for the sideward 
peak of IMF angular distribution, since it is not an easy task to analyze this 
phenomenon theoretically. 

In the first stage of the high-energy proton induced reaction, the incident 
proton interacts with target nucleons and produces several pions. In the inci- 
dent energy range around 10 GeV, the reaction mechanism of hadron-hadron 
collision evolves from single resonance baryon production to mutual excita- 
tion of incident baryons and/or string formation. For a reliable description 
of high-energy pA reactions, it is necessary to invoke these incident-energy 
dependent elementary cross sections parametrized and verified in pp and pA 
reactions. In Ref. [13], a hadron-string cascade model (JAM) is developed, 
and it is successfully appUed to hh and pA reactions as well as AA reactions. 
In the next stage of the reaction, scattered nucleons evolve in the mean field 
and then form an excited residual nuclei. Since the nucleon kinetic energy 
ranges from a few MeV to around 10 GeV and the fragment mass ranges from 
two (deuteron) to around 200, it is necessary to invoke appropriate mean 
field which takes account of the energy dependence as well as the density de- 
pendence and is verified in pA reactions. In Ref. [14], a parameterization of 
nuclear mean field is proposed, and this mean field is verified in pA reactions 
in a wide range within a framework of the Quantum Molecular Dynamics 
(QMD). In addition, it is already shown to well describe the incident energy 
dependence of the nucleon potential depth and the binding energy behavior 
from light to the heaviest nuclei. In the third stage of the reaction, the residual 
nucleus breaks up into fragments having various masses including IMFs, and 
the fragmentation mechanism is not understood well in the case we cannot 
assume that the system decays statistically from an equilibrated system. In 
heavy-ion reactions, there are some dynamical transport models which do not 
assume equilibrium [15-17] and well reproduce the IMF multiplicities. In this 
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case, the system expands to sub- saturation density due to the rapid collec- 
tive expansion from the initial hot and compressed nuclei. At sub-saturation 
densities, local density fluctuation is energetically favored, then fragments are 
easily formed in the dynamical stage. However, in the case of pA reactions, 
the flow effect is not large and the expansion is not enough to put the system 
into sub-saturation densities. At around the saturation density, density fluc- 
tuation is not energetically favored, but smoothed out during the evolution 
in a usual mean held type dynamics even if it is generated in the flrst stage 
of the reaction. Thus at present, there is no satisfactory microscopic model 
to describe the non-equilibrium fragmentation process in high energy proton 
induced reactions. 

In order to develop a model of non-equilibrium fragmentation in proton- 
induced reactions, we start from an equilibrium fragmentation model, percola- 
tion. It has been shown that the percolation model gives a good description of 
mass and momentum distribution of fragments in high-energy proton induced 
reactions, if we ignore the anisotropy of fragment emission in the residual 
nuclear rest frame. Since the initial nucleons are assumed to be distributed 
uniformly in a sphere in the equilibrium percolation model, the fragment an- 
gular distribution becomes automatically isotropic in the residual nuclear rest 
frame. Therefore, it is expected that the fragment angular distribution would 
be anisotropic if the initial nucleon distribution contains dynamical informa- 
tion of the earlier stage. We here formulate a new Non-Equilibrium Percolation 
(NEP) model which takes account of initial density and momentum fluctuation 
generated in the earher dynamics. In NEP, the nucleon phase space variables 
are taken from the results of the transport model calculation as the initial con- 
dition of percolation, instead of putting nucleons on sites. The bond breaking 
or connection probability is calculated by using the distance and momentum 
difference between the two nucleons under consideration, instead of giving a 
common bond breaking probability for the bonds connecting nearest neigh- 
bor sites. By doing these modiflcations, we can take into account the position 
and momentum fluctuations of nucleons inside fragment source and study the 
effects of non-equilibrium features of residual nucleus to the IMF formation 
processes. 

In this paper, we study the sideward enhanced IMF formation mechanism 
by using a combined framework of a transport model and a newly devel- 
oped Non-Equilibrium Percolation model. In the transport part, we use the 
cross sections and particle production algorithm developed in JAM [13], and 
the mean held verifled in QMD [14]. Hereafter we call this transport model, 
JAM/MF (JAM with mean held). The nucleon phase space information given 
in the transport model is used as the initial condition of NEP. We find that 
we can well reproduce the mass and energy distributions of IMFs in a unified 
way by using the above combined framework. In addition, if we assume that 
fragments are produced within a short time scale, we can also describe the 
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sideward enhancement of IMF emission. 

We find the foUowing reaction mechanism to enhance sideward emission of 
IMF: In the short time scale (around 20 fm/c), the dynamical effect in the first 
hadronic cascade stage still remains. Nucleons along the path of the leading 
incident proton are hcated-up and have large momenta in average. Since IMFs 
arc formed mainly in the cold region, the probability of IMF formation is 
small around the leading proton path. Then the IMF formation points are 
distributed non-spherically, predominantly distributed in a doughnut shape. 
After forming fragments in this doughnut shape, Coulomb repulsion between 
fragments pushes more strongly in the sideward direction in the rest frame of 
the residual nuclei. 

The outhne of this paper is as follows. In Sec. 2, we explain the transport model 
JAM/MF used in this work. In Sec. 3, a newly developed Non-Equilibrium 
Percolation (NEP) model is explained, and we discuss the validity of this 
model to show the good reproduction of fragment mass distribution. Next in 
Sec. 4, we discuss the mechanism of sideward enhancement of IMF emission 
through the analysis of the production position distribution of IMF, the IMF 
energy distribution and the IMF angular distribution. Finally in Sec. 5, we 
summarize this paper. 



2 Transport Model for High Energy Proton Induced Reactions 

In high-energy {Ep>j lOGeV) proton induced reactions, there arc mainly three 
stages — hadronic cascade, mean field evolution, and fragmentation. During 
these stages, nucleon energy decreases by about four order of magnitude from 
around 10 GeV to a few MeV. A large part of the proton incident energy 
(~ 10 GeV) is carried away by produced pions and emitted high momentum 
nucleons, and only a part of this energy remains as the excitation energy of the 
residual nucleus (several hundred MeV to a few GeV). This residual nuclear 
excitation energy is mainly exhausted to breakup the residual nucleus into 
fragments and nucleons, then the energy of finally produced fragments are 
around a few MeV per nucleon. 

For a reliable microscopic description of dynamical evolution of the system 
throughout these stages, it is necessary to apply reliable hadronic cross sections 
for hadronic cascade processes, and well verified mean field. In this work, we 
have developed a transport model, JAM/MF, by extending the hadron-string 
cascade model JAM [13] to include the mean field developed in Ref. [14]. It can 
be also considered as a version of QMD with hadronic cross sections developed 
in JAM. The results of JAM/MF are used as the initial configuration of the 
Non-Equilibrium Percolation described in the next section. 
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2.1 Hadronic cascade processes in JAM 

Here we describe the outline of the hadron-string cascade model, JAM. The 
detail of JAM is described in Ref. [13]. The main components of JAM are as 
follows. 

(1) Nuclear collision is assumed to be described by the sum of independent 
binary hh coUisions. Each hh collision is realized by the closest distance 
approach. In the original version of JAM, no mean field effect is included, 
therefore the trajectory of each hadron is straight in between two-body 
collisions, decays or absorptions. 

(2) The initial position of each nucleon is sampled by the parameterized 
distribution of nuclear density. Fermi motion of nucleons are assigned 
according to the local Fermi momentum. 

(3) All established hadronic states, including resonances with masses up to 
around 2 GeV, are explicitly included with explicit isospin states as well 
as their anti-particles. All of them can propagate in space-time. 

(4) The inelastic hh collisions produce resonances at low energies while at 
high energies {■\/s> 4GeV in BB collisions ^/s> 3GeV in MB coUisions 
and ^/s^2GeV in MM collisions) color strings are formed and they decay 
into hadrons according to the Lund string model [18]. Formation time is 
assigned to hadrons from string fragmentation. Formation point and time 
are determined by assuming yo-yo formation point. This choice gives the 
formation time of roughly 1 fm/c with string tension k — IGeV/fm. 

(5) Hadrons which have original constituent quarks can scatter with other 
hadrons assuming the additive quark cross section within a formation 
time. The importance of this quark(diquark)-hadron interaction for the 
description of baryon stopping at CERN/SPS energies was reported by 
Frankfurt group [19,20]. 

(6) At very high energies lOGeV), multiple minijet production is also 
included in the same way as the HIJING model [21] in which jet cross 
section and the number of jet are calculated using an cikonal formalism for 
perturbative QCD (pQCD). Hard parton-parton scatterings with initial 
and final state radiation are simulated using PYTHIA [18] program. 

(7) Pauli-blocking for the final nucleons in two-body collisions are also con- 
sidered. 

(8) We do not include any medium effect such as string fusion to rope [22,19], 
medium modified cross sections and in- medium mass shift. All results 
which will be presented in this paper are those obtained from the free 
cross sections and free masses as inputs. 

Hadron (proton, rr'^ ,tt~ , , K^) transverse mass rrit in proton induced re- 
actions (p(14.6GeV/c) + Be, Al, Cu, Au), light-heavy ion induced reactions 
(Si(14.6A GeV/c)+ Al, Cu, Au), and heavy-ion reactions (Au(11.6A GeV/c) 
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+ Au) has been analyzed in JAM, and it is found that JAM results show good 
agreement with the experimental data without any change of model parame- 
ters [13]. This means that momentum distributions of isolated single proton 
are correctly calculated by coUision cross sections of JAM. 

On the other hand, it is required to take account of later mean field evolution 
and percolation processes for the description of IMF production. 



2.2 Propagation in Quantum Molecular Dynamics 



In the framework of Quantum Molecular Dynamics (QMD), the classical tra- 
jectory of hadrons is determined by the stochastic hadronic collisions, de- 
scribed in the previous subsection, and the Newtonian equation of motion. 
The single particle wave function of each nucleon is represented by a Gaussian 
wave packet, having the phase space centroid parameter of (Rj, Pj) for the i-ih. 
nucleon. The total wave function is assumed to be a product wave function of 
nucleon Gaussian wave packets. Then wc can derive the Newtonian equations 
of motion on the basis of the time dependent variational principle [23]. The 
equations of motion become the canonical equations of the Gaussian centroid 
parameters as phase space variables, 

R.-^, P^--^- (1) 

The Hamiltonian H consists of the kinetic energy and the effective interaction 
energy. 



H = T + V, (2) 



^=E^^ = Ev^? + p?' (3) 

i i 

y — Vskyrme + ^Sym + Hlom + VcomI + ^Pauli ■ (4) 

In this interaction energy, the following terms are included: Skyrme type den- 
sity dependent interaction (Vskyrme); symmetry energy (T^ym), momentum de- 
pendent interaction taking account of the energy dependence of the nucleon 
potential (VMom), Coulomb potential between protons (Vcouiomb), and the Pauli 
potential (Vpauii)- The last term is added to take into account the Pauli prin- 
ciple approximately. The form of each term is shown in the Appendix. 

In this study wc use parameters shown in Table A.l. These parameters are de- 
termined to reproduce various kinds of systematic observables such as (l)the 
density dependence of the total kinetic energy of Fermi gas, (2)density depen- 
dence of the total energy of symmetric nuclear matter above normal nuclear 
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density, and (3)energy dependence of the real part of proton-nucleus optical 
potential, binding energy and radius of heavy nuclei {A > 3). All these observ- 
ables are calculated by metropolis sampling method under the condition that 
the temperature of nucleus is 2.5 MeV. Including the Pauli potential improves 
the reproductivity of the quasi-elastic part of the neutron energy spectra in 
p(256MeV,800MeV)+2°«Pb [14]. 

In the simulation of 12 GeV proton induced reaction, the ground state of 
target nucleus (Au) is calculated by frictional cooling and heating method to 
fit the binding energy of Au. 

In the total Hamiltonian H we have already introduced the relativistic form 
of kinetic energy expression and adopted the relativistic kinematics in the 
cascade processes in JAM. In order to make a full relativistic description, 
it is also necessary to formulate the interaction terms in a covariant way. A 
Lorentz-covariant extension of the QMD, debbed relativistic quantum molec- 
ular dynamics (RQMD), has been proposed in Ref. [24], which is based on the 
Poincare-invariant constrained Hamiltonian dynamics. Although the RQMD 
is a numerically feasible extension of QMD towards a fully covariant approach, 
it still costs too much computing time to apply to heavy systems. Therefore, 
we make an alternative extension of QMD [25]. We introduce Lorentz-scalar 
quantities into the arguments of the interactions. 



In these expressions, thin letters {Dij, Qij, Pij) show four vectors, and Pij 
denotes the two-body center-of-mass four momentum 



It is easy to show that the quantity D^^ is the two-body distance squared in 
the two-body CM system. By these procedures we can take into account the 
primary part of the relativistic dynamical effects approximately in our QMD. 



3 Fragment Formation 

There are many attempts to describe the IMF formation in multifragmenta- 
tion processes by using microscopic molecular dynamics (MD) simulations. 
Since nuclear fragmentation is a truly quantum-mechanical process and it is 
beyond the applicability limit of semiclassical MD theory, it is necessary to 
improve MD theory to include quantum mechanical effects [16,17,26-29]. In 



Dj = (R, - R,)2 ^ Dj = -Dl + (/^, ■ A,)V4 ' 

= (P^ - P^)' - = -Ql + {P^^ ■ Q^jf/Pi ■ 



(5) 
(6) 



Pij = Pi + Pj = {Ei + Ej, Pi + P,) . 



(7) 
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Ref. [26,29], the widths of nucleoli wave packets are treated as time-dependent 
dynamical variables. With this extension, we can find significant improve- 
ments in ground state properties such as binding energies, density profiles 
or a-clustering structure in light nuclei, nucleon emission or disintegration of 
clusters in medium energy heavy-ion collisions. Another type of proposed mod- 
ification is to include quantum fiuctuation or quantum branching [16,17,27]. In 
Rcfs. [16,27], the authors discussed the necessity of quantum fiuctuation in the 
equation of motion of MD in order to ensure that the wave packet distribution 
follows the quantum statistics at equilibrium. The effects of quantum fluctua- 
tion in fragmentation processes such as evaporation of particle and fission like 
fragmentation are discussed in Ref. [27]. The necessity of quantum branching 
is discussed in Rcfs. [17,28]. The authors introduce the quantum branching 
process to compensate the restriction on the one body nucleon distribution, 
and incorporate the diffusion and the deformation of wave packets accord- 
ing to the Vlasov equation [17]. By this modification, the nucleon motion in 
the mean field is improved, and the fragment mass distribution is successfully 
reproduced in heavy ion reactions. 

It may be possible to improve the description of IMF formation also in high- 
energy pA reactions by making above modifications. However, all of these 
modifications are formulated in non-relativistic ways, and the way to include 
quantum effects is under discussion. Thus in this work, we treat the fragmen- 
tation process in a more phenomenological manner. 

In a standard treatment of fragment formation, first we calculate the forma- 
tion of excited fragments in a transport model or in a phenomenological way, 
then provided that these fragments are thermalized they are forced to de-excite 
through particle emission and/or fission in a statistical model by using only a 
few quantities of fragments such as the mass, excitation energy and angular 
momentum. Since we would like to address the non-equilibrium properties of 
fragment formation in high-energy pA reactions, here we do not assume that 
nucleons are distributed uniformly but utilize information in the dynamical 
stage — nucleon spatial and momentum distribution inside the residual nu- 
clei generated in the transport model, JAM/MF. In order to activate these 
nucleon phase space degrees of freedom in fragmentation processes, the perco- 
lation model would be an appropriate starting point. In percolation models, 
nucleon positions are already considered, and we can easily assign the momen- 
tum to each nucleon. 

In this section, we explain our prescription to form residual nuclei from trans- 
port model results, and how to incorporate non-equilibrium nature in the 
percolation model. 
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3. 1 Residual nucleus formation 



After successive collisions and mean field evolution of nucleons, several parti- 
cles arc emitted and a residual nucleus is formed. Here, we take all phase space 
information of all particles consisting the system. Then we can obtain physi- 
cal property of residual nucleus such as charge, mass and excitation energy as 
follows. 

We assume an external nuclear mean field of nucleons around the residual 
nuclei having a depth of 50 MeV. Then we assume that particles are trapped 
and form a residual nucleus when the nucleon single particle energy is negative 
(for neutrons) or below the Coulomb barrier (for protons). The excitation 
energy E* of the residual nucleus is calculated in a similar way to that in the 
exciton model [30] . Here, we assume the nucleon Fermi energy in the residual 
nuclei £f is 40 MeV. 



^* = E4 + E4' (8) 

i=l 3=1 

Si ^Ep-Si , e]^ Ej - Ef . (9) 

Here Ei denotes the target nucleon single particle energy just before the first 
collision of that target nucleon, then the hole number (A^^) corresponds to the 
number of target nucleons which collide with the leading proton or secondary 
cascade particles. The particle energy Ej is the scattered nucleon single particle 
energy at the time of residual nucleus formation. The number of particles {Np) 
is calculated at the residual nucleus formation. This method is basically the 
same as that of Ref. [31]. 



3.2 Non-Equilibrium Percolation model 



In (bond) percolation models, nucleons in the residual nucleus are put on the 
lattice sites, and we cut the bond between the nearest neighbor sites accord- 
ing to the bond breaking probability given as a function of excitation energy. 
Then the connected nucleons are assumed to form fragments. By estimat- 
ing the excitation energy and the mass of the residual nucleus in a simple 
physical model, it has been demonstrated that the percolation model works 
very well to describe the mass and energy distribution of fragments produced 
from multifragmentation processes in relativistic energy ( > 10 GeV ) pA 
reactions [32-36]. 

In standard percolation models, there is one strong assumption which disables 
us to describe IMF anisotropy: the total excitation energy is assumed to be uni- 
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formly distributed over the entire excited system of the residual nucleus, and 
the nucleons are spatially distributed uniformly in a sphere. This is equivalent 
to assume that equilibrium is reached in the residual nucleus. This assumption 
is very convenient to discuss some formal aspects of the model, such as the 
critical bond breaking probabilities, critical exponents, and so on. In reality, 
the nucleon may not be uniformly distributed spatially and the excitation may 
depend on the position inside the residual nucleus, if the fragmentation time 
scale is short. Therefore, as a natural extension of the percolation model to 
include non-equilibrium effects, we start percolation from nucleon phase space 
distribution given in the dynamical transport model calculation. In this case, 
since nucleons are not necessarily on the lattice sites, the bond breaking prob- 
ability becomes a function of distance between the nucleon pair. In addition, 
we can consider the position dependence of the bond breaking probability. 

Here we propose a Non-Equilibrium Percolation (NEP) model without lattice 
combined with a transport model. We start from nucleon phase space distri- 
bution given in the microscopic transport simulation, and the bond breaking 
probability is parametrized as a function of nucleon pair distance and posi- 
tions. In the following, we explain the actual calculation procedure in NEP. 

(1) Dynamical evolution of the pA reaction is simulated by using the trans- 
port model described in the previous section, Sec. 2, until t = tsy^. 

(2) Phase space information of nucleons inside the residual nucleus, and the 
bulk quantity of the residual nucleus such as charge 

^rest mass Aj-gg and 

excitation energy E* are obtained according to the prescription described 
in the previous subsection, Subsec. 3.1. 

(3) All nucleons inside the residual nucleus are connected by bonds. 

(4) Break the bond when the bond length L5 is large enough. When the total 
excitation energy exceeds a given value, E* > Eo{Ares) {Eq = SciAres 
MeV), the bond is broken if > Rcut- When the total excitation energy 
is small, E* < Eo{Ares): bonds are broken if > C2-Rcut- This total 
excitation energy dependence is necessary to keep the residual nucleus 
stable at low excitation. 

(5) The bond between the i-th and j-th nucleons is broken with a probabil- 
ity Pb{i,j)- The bond breaking probability Pb{i,j) depends on the local 
excitation of the residual nucleus, and it is parametrized as follows. 



Here, E*(Rij) is the local excitation energy at the center-of-mass position 
(Rij = (Ri-|-Rj)/2) of the nucleon pair. We calculate E*{Iiij) as follows. 



Pb{i,j) = 1 - exp [-iyE*{Rij)] 



(10) 




(11) 
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* ^ / tfc - e (tfc > e) 
' 10 {h<e) 



(12) 



where tk is the kinetic energy of k-th nucleon in the CM frame of the 
residual nucleus and is the position of k-th nucleon. 
(6) Fragments are identified as the bond connected nucleons at this stage. 
The fragment momentum Pf and position Rj are determined as the CM 
momentum and position of nucleons belonging to this fragment. 

The model parameters, Rcut,^,o,,e,Ci,C2, are determined to fit the mass dis- 
tribution of fragments produced in p(11.5GeV)-|-Au reaction. The values of 
NEP parameters are shown in Table 1. 

Table 1 

Parameters of NEP 
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Fig. 1. Calculated fragment mass distributions in p(12 GeV)+Au reaction. Dashed 
histograms and solid histograms show calculated results of JAM/MF and JAM/MF 
followed by NEP, respectively. In the left, middle and right panels, the dependence 
on the the switching time {tsw =15, 20, and 30 fm/c, respectively) is shown. Solid 
points show the experimental data of p(11.5GeV)-|-Au reaction [37]. In the transport 
model, simulation calculations are made in the impact parameter range of 6 < 7.5 
fm. 

In Fig. 1, we show the calculated fragment mass distributions in p(12 GeV)+Au 
reaction in comparison with the experimental data [37]. The mass of residual 
nucleus (dashed histograms) is distributed in the range 160<, Ares^ 197. This 
means that, in the transport model simulation, the number of emitted nucleons 
is less than 40. Then, the decay of the residual nucleus to various mass frag- 
ments is described in NEP as shown by the solid histograms. One can see that 
the characteristic U-shape curve of the experimental mass distribution is well 
reproduced in our calculation, and this behavior does not strongly depends on 
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the switching time. When we analyze the impact parameter dependence of the 
mass distribution, we find that this U-shape distribution is a consequence of 
the integration over the impact parameter. At central collisions, the residual 
nucleus generally has smaller mass and larger excitation energy, and it decays 
into smaller fragments, in average. At peripheral collisions, the residual nu- 
cleus is not highly excited, and emits several nucleons even in the percolation 
stage. 

It should be noted that the above global feature of the mass distribution can be 
also reproduced in a standard percolation model. In Fig. 2, we show the model 
dependence of the mass distribution. Here we compare the combination of the 
models, JAM with Statistical Decay Model (SDM), Percolation model, JAM 
with Percolation model, and JAM/MF with NEP. In SDM, the evaporation of 
light fragments {n,p, d, t, "^He, a) and the fission of heavy nuclei are included, 
but these are not enough to form IMFs abundantly. On the other hand, in the 
percolation model of Ref. [35] which fit the mass distribution of p+Ag reaction 
at 300 GeV, the global mass distribution is reproduced without any change 
of the model parameters, and it is possible to get a better fit by using target 
dependent parameters. 

This comparison of model calculations suggests the necessity to invoke some 
kind of multifragmentation mechanism such as the percolation in addition to 
the evaporation and fission. If we adopt the percolation model in fragmenta- 
tion stage, we can qualitatively describe the U-shape mass distribution, which 
we cannot explain in statistical decay models with evaporation and fission. 
Furthermore, in order to reproduce the whole range of the mass spectrum, it 
is necessary to stabilize the nucleus at low excitation. In the original percola- 
tion model [35] (upper right panel), nucleon positions are kept to be on sites, 
and the bond breaking probability becomes smaller in peripheral collisions, 
then larger mass nuclei can survive. With QMD transport (lower right panel), 
heavy residual nuclei can survive percolation at low excitation, because the 
mean field keeps the nucleus stable. Without mean field (lower left panel), on 
the other hand, nucleons with moderate kinetic energy can easily escape from 
residual nuclei, then we cannot form stable heavy nuclei. 

This comparison tells us the merit to use a combined framework of JAM/MF 
with percolation for the analysis of the fragmentation processes. However, it 
is possible to reproduce the mass spectra without invoking spatial and/or 
momentum fiuctuation generated in the dynamical stage by modifying pa- 
rameters. Therefore, for the investigation of non-equilibrium nature in the 
fragmentation, it is necessary to analyze the energy and angular distribution 
of fragments. 
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Fig. 2. Model comparison of fragment mass distribution in p(12 GeV)+Au reaction. 
Four model results — JAM+SDM (upper left), Percolation model with default 
parameters (upper right) [35], JAM+Percolation (bottom left), JAM/MF + NEP 
(bottom right) — are shown in histograms and compared with the experimental 
data [37]. As for the explanation of these models, see the text. Dashed histograms 
in the lower panels show the calculated results multiplied by 1.3 and 1/1.3, showing 
that the combined analyses of a transport model with percolation reproduce the 
data within around 30 % accuracy. 



4 Results 
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4 . 1 IMF formation points 



Now we analyze the non-equilibrium properties of IMF production in the 
framework explained in the previous sections. In Fig.3, IMF(6 < Z < 20) 
formation position distributions in xy-plain calculated in the transport model 
calculation followed by NEP are shown. In the transport model simulations, 
5000 events are generated, and the results of NEP from these events are over- 
layed in the figure. In central events (6 < 3 fm, upper panels), we can see that 
IMFs are produced in the doughnut shaped region. Along the leading proton 
path, a large part of nucleons collide with the leading proton or secondary 
cascade particles. Since they have large kinetic energies, they are not judged 
to belong to the residual nuclei. Even if they stay in the residual nuclei, the 
local bond breaking probability becomes large. As a result, the IMF forma- 
tion is suppressed along the leading proton path. In the surrounding doughnut 
shaped region, the number of collision is small, then the excitation energy is 
low. 



In the lower panels of Fig.3, IMF formation points in peripheral collisions (3 
fm < 6 < 7.5 fm) are shown. We can see that IMFs are produced from a new 
moon shaped region. In peripheral collisions, highly excited region is shifted 
to outer side. Then peripheral excited region shrinks to narrower at time step 
30fm/c. 



IMF production position distribution 
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Fig. 3. Calculated IMF {6 < Z < 20) position distributions in xy-plain. Solid points 

display the IMF formation points in NEP following the transport model evolution 
until the switching time of 15, 20 and 30 fm/c. 5000 events of the transport model 
events are overlayed. In upper and lower panels, central (6 < 3 fm) and peripheral 
(3 < 6 < 7.5 fm) events are shown, respectively. 
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In Fig. 4, we show calculated IMF formation points in ,T2;-plain, restricted to 
\y\ < 2.0 fm. We can see that IMF formation point distributions clearly have 
a hole along the leading proton path in central collisions (upper panels) and 
at early times in peripheral collisions (lower panels). 

In both of the cases, IMFs are formed mainly on the surface of the residual 
nuclei, and the formation is suppressed in the central part. This feature can 
be understood as follows. When the excitation energy is small, nucleons are 
hkely to be connected with each other. In the central part of the residual 
nucleus, IMF can grow easily to heavier nucleus by being connected with cold 
nucleons in the surrounding region. This mechanism may emerge also in the 
equilibrium percolation model. However, this will not generate any anisotropy 
in IMF angular distribution. 

IMF production position distribution 



15 



t=15fm/c 

b < = 3 fm 


. t=20 fm/c 

■X- 

b < = 3 f m 


, t=30 fm/c 

..-'^ : -:■ . 
••.fcv-.." 

b < = 3 fm 


t=15fm/c 

• . 

3 fm < = b <= 7.5 fm 


t=20 fm/c 

3 fm < = b < = 7.5 fm 


t=30 fm/c 

•f • • 

3 fm < = b < = 7.5 fm 



X 

Fig. 4. Calculated IMF {6 < Z < 20) position distributions in xz-plain, restricted 
to \y\ < 2.0 fm. Solid points display the IMF formation points in NEP following 
the transport model evolution until the switching time of 15, 20 and 30 fm/c. 5000 

events of the transport model events are ovcrlaycd. In upper and lower panels, 
central (6 < 3 fm) and peripheral (3 < 6 < 7.5 fm) events are shown, respectively. 

4-2 IMF energy distribution 

After the fragments being simultaneously produced through the multifragmen- 
tation of the residual nucleus, the fragments are dispersed by the Coulomb 
repulsion. In the uniform Coulomb expansion of fragments, the total frag- 
ment kinetic energy is calculated as the sum of the thermal and the Coulomb 
energies. In order to describe the Coulomb expansion of fragments, here we 
consider the following classical Hamiltonian H, 
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i=l i<j 



R 



(13) 



where m/j, R/j, and P/j are the mass, position, charge, and momentum of 
the i-th fragment, respectively. Prom this Hamiltonian Eq.(13), the Newtonian 
equations of motion arc solved numerically. The initial values of the R/j and 
P/j are given in NEP as discussed in Sec. 3. The time evolution of the Coulomb 
expansion is solved until the total kinetic energy of each fragment ceases to 
show appreciable changes. Then we obtain kinetic energies of fragments in the 
final state. Here, we adopt the time step of Ai=0.5 fm/c and integration until 
500 time steps are carried out. 

Energy Distribution of Fragment 
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Fig. 5. Calculated energy distributions of O, F and Ne isotopes produced in p{Vl 
GeV)+Au reaction in comparison with the experimental data. The vertical axis 
shows the cross section and the horizontal axis shows the kinetic energy of the 
fragment. Left, middle and right panels show the energy distributions of O, F and 
Ne, respectively. Upper, middle and lower panels show the energy distributions 
emitted at angles of 30, 70 and 110 degrees from incident proton direction in the 
laboratory frame. Solid and dashed histograms show the calculated results with and 
without the Coulomb expansion process. The results of switching time tgw = 20 fm/c 
are shown. Solid points are the experimental data [38] 

In Fig. 5, we compare the calculated IMF energy distributions produced in 
p(12 GeV)+Au reaction with the experimental data [38]. Prom these results 
we can see large effects of Coulomb expansion in the energy distribution of 
IMF. If we take into account only the kinetic energy of fragments calculated 
in the transport model followed by NEP, the energy peak around 50 MeV is 
never reproduced. On the other hand, we find that the calculated results after 
Coulomb expansion reproduce the qualitative behavior of the experimental 
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data well; peaked around 40-60 MeV, enhancement at sideward angles, and 
the absolute values of the double differential cross sections. These facts show 
us that the Coulomb expansion process plays a very important role in the 
fragmentation phenomenon. 

4-3 Analysis of IMF angular distribution 

In Fig. 6, we show the calculated switching time dependence of the IMF angu- 
lar distribution. In comparing the results with (left panel) and without (right 
panel) Coulomb expansion, the Coulomb expansion effects are found to modify 
the angular distribution largely from forward peaked to sideward peaked. This 
mechanism is connected with the IMF formation point in the residual nucleus. 
As shown at Fig. 3 and Fig. 4, IMFs are produced dominantly in a doughnut 
shaped region in violent central collisions. Then the Coulomb force among 
fragments pushes IMFs sideways of the doughnut region even if IMFs are for- 
ward peaked before the Coulomb expansion. In peripheral coUisions, although 
IMFs are formed in a new moon shaped regions at early times, the rapid en- 
ergy transport washes out the anisotropy of the IMF production points. Then 
if the impact parameters are integrated over, the angular distribution of IMFs 
becomes forward peaked at around the switching time tsw = 30 fm/c even after 
the Coulomb expansion as shown in the left panel of Fig. 6. This observation 
is consistent with the behavior expected in a equilibrium statistical picture of 
fragmentation. Therefore, our present study suggests the possibility that IMF 
angular distribution would have sideward peak at around 90 degrees in the 
rest frame of the residual nuclei if IMFs can be formed in a very short time 
scale around 20 fm/c. 



5 Sumrriciry 

In this paper, we have analyzed the mechanism of formation and sideward 
enhancement of IMFs produced in the 12 GeV proton induced reaction on 
Au target by using a combined framework of a transport model (JAM/MF) 
followed by a newly developed Non-Equilibrated Percolation (NEP) model. 
NEP can take into account non-uniform spatial distribution of nucleons and 
the local excitation in the residual nucleus, then it would be appropriate to 
analyze non-equilibrium natures of IMF formation. The present combined 
model has been verified in various ways. In the hadronic cascade stage (JAM), 
momentum distribution of nucleons has been verified by comparing hadron 
(proton, 7r+,7r~,i^"'",i^'~) transverse mass spectra [13]. The effective interaction 
in QMD has been verified through the study of matter and nuclear proper- 
ties and the quasi-elastic neutron energy spectra in pA reactions at 256 MeV 
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Angular Distribution of IMF 
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Fig. 6. Calculated switching time dependence of the IMF angular distribution. 
Dotted, dash-dotted, dashed and solid histograms display the results at switching 
time of 15, 20, 22.5 and 30 fm/c, multiplied by 0.4, 0.6, 0.75 and 1.0, respectively. 
In the left (right) panel, the results with (without) Coulomb expansion is shown. 

and 800 MeV [14]. In the fragmentation (NEP), although there are still sev- 
eral model parameters, they arc determined to fit the fragment mass spectra 
in high-energy pA reactions. By using the same parameter, we applied our 
model to the fragment energy distributions, and the experimental data are 
well reproduced. 

By using this model we have shown that IMFs are produced in a doughnut 
shaped region and a new moon shaped one inside the residual nucleus in cen- 
tral and peripheral collisions, respectively. We suggest that combined effects 
of such a geometry of IMF formation point distribution and the Coulomb 
expansion between fragments would lead to a sideward peaked angular dis- 
tribution of IMFs provided that the IMFs are formed in a very short time 
scale of around 20 fm/c. In principle, it would be possible to estimate this 
time scale experimentally by measuring the particle correlations in momen- 
tum multi-dimensions. Another remaining interesting problem is to study the 
mechanism of the incident energy dependence of the IMF angular distribution, 
which evolves from forward peaked, sideward peaked, to backward peaked as 
the incident energy grows. Our approach will be also useful in analyzing those 
problems. 
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A Quantum moIeculEir dynamics ( QMD ) 



In QMD method, the single particle wave function of the i-th nucleon is rep- 
resented in a Gaussian wave packet, having the density distribution. 



A(r) = 



1 (r-Rj)^ 

21, 



(27rL) 



(A.l) 



where L denotes the width parameter for the Gaussian wave packets and Rj 
are the centers of position of the i-th nucleon. 

Each term in the interaction energy in Eq. (4) depends on the center of position 
R of the nucleon and the center of momentum P of the nucleon as follows. 
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(A.2) 

(A.3) 
(A.4) 

(A.5) 
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Here, erf denotes the error function, and q is 1 for protons and for neutrons. 
We only take into account the effective interactions which affect nucleons. 
< Pi > is an overlap of density with other nucleons defined as 



< pi>= J2 = j d^Pii^)pji^) 
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1 



(A.7) 



< Pi >' means that L of < > is modified to be L' = i(l + t)^L. The 
parameters are listed in Table A.l. 

Table A.l 

QMD parameters 



Cp,,ai 110 [Mi^] 


25.0 [M(-V] 


qo IMA [fm] 


yii^ -258.54 [MeV] 


Po 120 [MeV/c] 


/xi 2.35 [fm-i] 


L 1.75 [fm2] 


Fex^ 375.6 [MeV] 


a -127.86 [MeV] 


112 0.4 [fm-i] 


P 204.28 [MeV] 


/9o 0.168 [fm-3] 


1 
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